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A generalized mathematical model is investigated of sloshing dynamics for dewar containers , partially filled 
with a liquid of cryogenic superfluid helium II, driven by both gravity gradient and jitter accelerations 
applicable to two types of scientific spacecrafts , which are eligible to carry out spinning motion and/or slew 
motion to perform scientific observations during normal spacecraft operation. Two examples are given for 
the Gravity Probe- B (GP-B) with spinning motion , and the Advanced X-Ray Astrophysics 
Facility-Spectroscopy ( AXAF-S) with slew motion , which are responsible for the sloshing dynamics. Explicit 
mathematical expressions for the modelling of sloshing dynamics to cover these forces acting on the spacecraft 
fluid systems are derived. The numerical computation of sloshing dynamics will be based on the noninertial 
frame spacecraft bound coordinate , and we will solve the time -dependent three-dimensional formulations of 
partial differential equations subject to initial and boundary conditions. Explicit mathematical expressions of 
boundary conditions to cover capillary force effects on the liquid-vapor interface in microgravity environments 
are also derived. Results of the simulations of the mathematical model are illustrated. 

Keywords: cryogenic helium, slosh dynamics, microgravity liquid-vapor interface, capillary effect 


1. Introduction 

In order to carry out scientific experiments, some 
experimental spacecraft use cryogenic cooling for 
observation instrumentation and telescopes and super- 
conducting sensors for gyro read out. They also maintain 
very low temperatures, near absolute zero, for mech- 
anical stability. The approaches to both cooling and 
control involve the use of superfluid liquid helium II. In 
this study, mathematical models of sloshing dynamics 
associated with spinning and/or slew motions are 
investigated. To cover the spacecraft spinning and/or 
slew motions, the Gravity Probe-B (GP-B) ( Figure 1) and 
the Advanced X-Ray Astrophysics Facility-Spectro- 
scopy (AXAF-S) (Figure 2) spacecrafts have been chosen 
as examples in this study. Both the GP-B and the 
AXAF-S spacecrafts adopt the cooling and boil-off from 
a cryogenic liquid helium dewar as a cryogen and 
propellant to maintain the cooling of instrumentations, 
attitude control, and drag-free operation of the 
spacecraft. The potential problems for cryogenic liquid 
in the dewar container could be due to asymmetry in the 
static liquid helium distribution and to perturbations in 
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the liquid-vapor interface caused by slosh wave 
excitation driven by pointing control, machinery 
vibration, etc. 

For the cases of both the GP-B and the AXAF-S 
spacecrafts, cryogenic liquid helium II, at a temperature 
of 1.3°K, is used as the propellant. Due to its superfluid 
behavior, there is no temperature gradient in the liquid 
helium. In the absence of temperature gradients along 
the surface, which drive Marangoni convection, the 
equilibrium shape of the free surface is governed by a 
balance of capillary, centrifugal, and gravitational forces. 
Determination of liquid-vapor interface profiles based 
on computational experiments can uncover details of the 
flow that cannot be visualized easily or measured 
experimentally in a microgravity environment. 

Instabilities of the liquid-vapor interface can be 
induced by the presence of longitudinal and lateral 
accelerations. Slosh waves are thus excited that produce 
high and low frequency oscillations in the liquid 
propellant. A recent study 1 suggests that the high 
frequency accelerations may be unimportant in compari- 
son with the residual motions caused by low frequency 
accelerations. 

Time-dependent dynamical behavior of partially filled 
rotating dewars in reduced gravity environments was 
simulated by numerically solving the Navier-Stokes 
equations subject to initial and boundary conditions. 2-5 
At the interface between the liquid and the gas, both the 
kinematic boundary condition and the conditions for 
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Figure 1. Gravity Probe-B (GP-B) spacecraft on-orbit configuration. 


Z 



Figure 2. AXAF-S spacecraft coordinate systems with azimuth 
angle t^ E from the spacecraft mass center to the center of the Earth. 
Coordinate (x", t) for slew motion and coordinate (x, y, z) for 
fluid mechanics computations. 


tangential and normal stresses at the interface were 
applied. 6-10 The initial conditions were adopted from the 
steady-state formulations developed by Hung et al. 11-14 
Some of the steady-state formulations of interface shapes 
were compared with available experiments carried out 


by Leslie 15 in a free-falling aircraft (KC-135). The 
experiments carried out by Mason et al. 16 showed that 
classical fluid mechanics theory is applicable for 
cryogenic liquid helium in large containers. 

In the spacecraft orbit around the Earth, the direction 
of the azimuthal angle of the Earth toward the location 
of the spacecraft mass center varies from 0° along the 
rolling axis of the spacecraft to various directions; 
therefore three-dimensional calculations are a must. 

As the spacecraft moves along the orbit, any fluid 
capable of motion relative to the spacecraft is subject to 
the acceleration that arises from the gravity gradients of 
the Earth. 17-19 The interaction between the mass of fluid 
and the spacecraft mass due to gravity gradient 
accelerations 18 is capable of exciting slosh waves and 
disturbing the fluid system, thus inducing fluctuations in 
the viscous stresses and their moments with a related 
effect on the containers of the spacecraft. Furthermore, 
sources of residual acceleration of gravity jitter, including 
atmospheric drag on the spacecraft, background gravity, 
spacecraft attitude motions arising from machinery 
vibrations, spacecraft slew motion, thruster firings, crew 
motion, etc., are also capable of exciting slosh waves in 
the fluid containers. 

Considering the dynamical behavior of helium II in a 
large rotating cylinder, a mixture of the superfluid and 
the normal fluid without separation of the two fluids is 
accounted for in the model computation. 20 21 The 
density concentration of superfluid is a function of 
temperature, which is also true for the surface tension 
and viscous coefficient for helium II. 22,23 In this study, 
the theory of viscous Newtonian fluids is employed with 
modification of transport coefficients adjusted to account 
for normal and superfluid density concentrations which 
are functions of temperature. 22,24-26 
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2. Functions of scientific observations and 
spacecraft motions 

Both the GP-B and the AXAF-S spacecrafts are 
sun-synchronous Earth satellites orbiting at 650 km 
altitude directly over the poles. The functions of scientific 
observation for both the GP-B and the AXAF-S 
spacecrafts and their motions are as follows. 

2.1 The GP-B spacecraft 

The GP-B spacecraft ( Figure 1) is a relativity 
gyroscope experiment to test two extraordinary, 
universified predictions of Albert Einstein’s general 
theory of relativity. 27,28 By using gyroscopes (those with 
electrically supported spheres, spinning in a vacuum, and 
others that utilize the spins of atomic nuclei, circulating 
sound waves, and even circulating laser beams), the 
GP-B measures two distinct space-time processes, frame 
dragging and the geodesic effect, which gradually change 
its directions of spin. In these gyroscopes, the underlying 
principle is that rotating systems, free from disturbing 
forces, should remain pointing in the same direction in 
space. 

To comprise these functions of scientific observation, 
the GP-B stores its gyroscopes, telescope, probe mass, 
and other equipment in the center core of the dewar 
probe surrounded by cryogenic helium II liquid. The 
dewar container of the GP-B is spinning with a rotation 
rate of 0.1 rpm during normal operation. Because the 
telescope is constantly oriented on Rigel, gyro spin 
directions also point toward Rigel during the spacecraft 
motion around the polar orbit. There is no slew motion 
involved in the GP-B with its rotation. 

2.2 The AXAF-S spacecraft 

The AXAF-S spacecraft {Figure 2) and its sister 
spacecraft AXAF-I (I for imaging) are two spacecrafts 
restructured from the original AXAF design to carry out 
astrophysical observations. Equipped with the (micro- 
calorimeter) X-Ray Spectrometer (XRS), the AXAF-S 
provides high throughput, high resolution, nondispersive 
spectroscopy at the high AXAF X-ray energies — includ- 
ing the astrophysically important iron-K spectral region 
(above 6.4 keV) — and also permits some spatially 
resolved high resolution spectroscopy. AXAF-S com- 
prises a foil mirror (or possibly, a replication optic) 
telescope (4.7 m focal length), with XRS in the focal 
plane. With the baseline optical system, the AXAF-S 
provides important, unique capabilities for high through- 
put, high resolution (above 1 keV) spectroscopy of 
extended and point sources, and for some spatially 
resolved high resolution spectroscopy. 

Because of the unique capabilities for high resolution 
spectroscopy of point and extended sources, AXAF-S is 
eligible to carry out the observations of distance scale 
stars, active galactic nuclei, clusters of galaxies, 
supernova remnants, X-ray binaries, etc. In this study, 
spacecraft slew motion without spinning with a rotating 
axis is required for the AXAF-S to perform its scientific 
mission. 


3. Basic characteristics of gravity gradient and 
gravity jitter accelerations 

Any fluid element inside the on-orbit spacecraft fluid 
system is subject to the acceleration that arises from the 
gravity gradient of the Earth. 1 7-1 9,29 “ 33 Once the 
spacecraft orbit is fixed, the orbit period is determined 
and the basic structure of the gravity gradient 
acceleration can also be calculated. However, the gravity 
gradient acceleration acting on each fluid element inside 
an on-orbit spacecraft fluid system is different since it 
depends upon the distance of the location of the fluid 
element to the mass center of the spacecraft and its 
direction toward the location of the center of the Earth. 
This acceleration can only be calculated based on a 
noninertial frame of spacecraft-bound coordinates. Thus, 
the coordinate system shall be transformed from 
ordinary inertial frame coordinates to noninertial 
coordinates. 

3.1 Orbit motion of spacecraft 

Let us consider the cases of the GP-B and the AXAF-S 
spacecrafts, which are the Earth satellites orbiting at 
650 km altitude directly over the poles. The orbit period, 
r 0 can be computed from the following expression: 


where R E denotes the radius of the Earth (6373 km), R c 
is the radius of the circular orbit ( = R E -b h = 7023 km), 
h is the orbit altitude (650 km), and g 0 is Earth’s gravity 
acceleration (9.81 m/s 2 ). For the case of both GP-B 
and AXAF-S spacecrafts, the orbit period is t 0 = 
97.6 minutes and the orbit rate is n = 2k/t q — 1.07 x 
10“ 3 rad/s. 

As the spacecraft is orbiting around the Earth, the 
azimuthal angle of the Earth, \f/ E toward the location of 
the spacecraft mass center varies with respect to time. At 
time t = 0, the rolling axis of the spacecraft is aligned 
with the radial direction of the Earth’s center to the 
spacecraft mass center. It is assumed that the spacecraft 
rolling axis is linearly turning from 0° to 360° in the orbit 
period, t c , when the spacecraft is orbiting around the 
Earth. This is particularly true for the case of the GP-B 
spacecraft. Without the spacecraft slew motion, the 
azimuth angle (\(/ E ) can be defined as 

4> E„ = t ( 2 ) 

where t is the time measured from the instant when the 
direction of the spacecraft rolling axis is aligned with the 
radial direction of the spacecraft mass center to the center 
of the Earth. 

3.2 Slew motion of spacecraft 

In order to carry out wide-range observations, some 
scientific spacecrafts require slew motion with respect to 
their mass centers. This is particularly true for the case 
of the AXAF-S spacecraft. For the case of the spacecraft 
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slew motion, the azimuthal angle, shown in equation (2), 
shall be modified through the coordinate transformation 
of slew motion when the spacecraft is orbiting around 
the Earth. 

Let us assume that the slew motion starts with the 
center located at the mass center of the spacecraft. Let 
us choose cartesian coordinates (x", /', z") with the 
z"-axis along the axis of the dewar container ( Figure 7). 

At time t = 0, the radial vector r c from the center of the 
spacecraft to the center of the Earth lies on 
the x"-z" plane of the cartesian coordinate chosen 
{Figure 2). The azimuth angle i j/ E is defined as the angle 
between the radial vector r c and the z"-axis. The 
rotation matrices for spinning and/or slew motions along 
the x"-, /'-, and z"-axes can be expressed as 


"l 

0 

0 


COS C O y t 

0 

— sin (o v i 


COS (O z t 

sin oj z t 

o’ 

0 

cos co x t 

sin o) x t 

, 

0 

1 

0 

, 

— sin o) z t 

COS OJ z 

0 

0 

— sin co x t 

COS (O x 


sin co y t 

0 

COS t Oyt 


0 

0 

1 


respectively. Here, o) x , co y , and o) z denote angular velocity of slew and/or spinning motions along the x"-, 
/'-, and z"-axes, respectively. The radial vector r c in cartesian coordinates without slew and spinning 
motion is ( Figure 2) 

/ co = [sin i/^ Eo , 0, -cos <Ae J (3) 

For spinning motion along the z"-axis, slew motion along the /'- axis, and slew motion along the 
x"-axis, the radial vector r c is 



"l 

0 

1 


COS COyt 

0 

— sin ojyt 


COS 0) z t 

sin co z t 

0“ 


sin iA Eo 

f = 

' c - z,y, x 

0 

COS iO x t 

sin o> x t 


0 

1 

0 


— sin a ) z t 

cos c o z 

0 


0 


0 

— sin t o x t 

cos co x t 


sin u) y t 

0 

COS COyt 


0 

0 

1_ 


— cos lp Eo 


Mathematical formulations for the dynamics of 
spacecraft slew motion derived are quite general. In order 
to show a specific example in numerical simulation, 
spacecraft slew motion that operates at 90° in 600 s along 
the /'-axis is considered in this study. 

In addition to the modification of the azimuthal angle 
introduced by the spacecraft slew motion through the 
formulation of coordinate transformation shown in 
equations (3) and (4), accelerations are also induced on 
the fluid mass in the dewar container. Accelerations 
acting on the fluid particle in the dewar induced by the 
slew motion of the spacecraft with the coordinates fixed 
at the spacecraft center of the mass is as follows (Figure 
2 ): 

ft p = cox((oxR p ) + axR p + 2o>xv (5) 

where ft p denotes the acceleraton vector of the fluid 
particle in the dewar container relative the body 
frame of the spacecraft, co is the angular velocity of the 


spacecraft body frame; a is the angular acceleration of 
the spacecraft body frame, and v is the velocity of the 
fluid particle relative to the spacecraft body frame. The 
first, second, and third terms in the right-hand side of 
equation (5) denote centrifugal, angular acceleration, and 
Coriolis accelerations, respectively. 

As indicated earlier, let us assume that the slew motion 
starts with the center located at the spacecraft mass 
center. Cartesian coordinates (x", /', z") are chosen with 
the origin located at the spacecraft mass center. Let us 
also assume that the x"-z" plane intersects the » enter of 
the Earth and the spacecraft mass center. In other words, 
the azimuthal angle of the Earth toward the spacecraft 
mass center lies in the x"-z" plane. Slew motion is along 
both the x"- and /'-coordinates. Thus, o s = (u) sx , c o sy , 0) 
and a s = (a sx , a sy , 0), R p due to slew motion becomes 


X" 


a> s /cu SJC Ry R x co %y j *+■ ot sy R z + 2o) sy v z 


= 

— co sx (co sx Ry R x co sy ) ol sx R z 2oj sx v z 

R,_ 

slew 

_-R z (Wsx + ojfy) 4- (cc sx R y - oc sy R x ) 4- 2{(o sx v y - o> 9y v x )_ 
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3.3 Coupling of the accelerations due to spinning and 
slew motions of the spacecraft 

As indicated in Section 2 of this study, in specific 
missions of scientific observation, the dewar container of 
the GP-B is spinning with a certain rotating rate without 
slew motion during normal operation, while the AXAF-S 
requires slew motion for pointing control for observation 
of extended sources of astronomical objects without 
spinning. For some particular reasons required in other 
spacecraft, both spinning and slew motions may be 
needed simultaneously. In this case, the following 
expression with coupling for the accelerations of spinning 
and slew motions of the spacecraft holds: 


ft 


p, slew and spinning 


R* 

Ry 

p 

L 2 Jslew and spinning 

^syi^sx Ry ^syR-x) "F yRz "F 2Ci) s yV z 

t^sxi^sxRy xRz 2(D S x^z 

L Rzi^sx "F ^sy) "F ( a s X Ry &syRx) "F 2(CQ sx V y ^ S y^x) 


+ 


{o) z R x o} &x R 2 )o) z oj z Ry 2co z Vy 
(co sy R z — a) z R y )aj z + cb z R x + 2 <o 2 v x 
(^sx^jc "F ^sy^y)^z 


spinning and coupling 


where co z and d> x denote angular velocity and angular 
acceleration, respectively, of the spacecraft spinning 
motion along the z-axis. 

Using ( R x , R y ) = (r cos 0, r sin 0) and (r x , t? y ) = (u r • 
cos 0 — u 0 sin 0, u r sin 6 + u 6 cos 6) for spinning motion 
only, equation (7) converted from cartesian to cylindrical 
coordinates becomes 


( 7 ) 


ft 

^p, spinning 


and 


R* 


— r cos 6a>l — r sin 0d> z — 2 (u r sin 0 + u e cos 0) w z 


Ry 
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spinning 
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spinning 


R r 


R x cos 6 + R y sin 6 


— roj>\ — 2 u B co z 

Re 

= 

— K x sin 9 + R cos 9 

= 

rd> z + 2u r co z 

R z 

[spinning 


spinning 

0 


( 8 ) 


( 9 ) 


3.4 Gravity gradient acceleration 

The gravity gradient acceleration acting on the fluid 
mass of spacecraft is 

-" 2 [3(r c -d)r c - d] (10) 

where a gg denotes the gravity gradient acceleration 
vector, d is a vector (not a unit vector) from the fluid 
element to the spacecraft geometric center, f c is a unit 
vector from the spacecraft geometric center to the center 
of the Earth, and n is the orbit rate (Figure 2). 

It is assumed that the gravity gradient exerted on the 
geometrical center of the spacecraft orbiting around the 
Earth on its specified orbit is zero. In other words, all 
the gravity acceleration exerted on the spacecraft is 
nothing but the gravity gradient acceleration, which is 
defined in equation (10). 


For the convenience of mathematical calculations, let 
us describe all the parameters involved in equation (10) in 
terms of cartesian coordinates. In order to match with 
the computer simulation, mathematical derivations are 
considered in the first quadrant. Figure 2 illustrates the 
geometrical relationship of the parameters shown in 
equation (10). 

Let us consider the fluid element of interest, m, located 
at (r, 0, z) in cylindrical coordinates and at (x, y, z) in 
cartesian coordinates. The spacecraft mass center is 
located at z = L c . Assume that vector f c lies in the x-z 
plane of the cartesian coordinate system. Vector d in 
(x, y, z) coordinate becomes 

d = [ — r cos 0, — r sin 0, — (z — L c )] (11) 
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Substituting equations (4) with cd x = oj y = 0 and 
equation (11) into equation (10), the noninertial frame 
expression of the gravity gradient acceleration with 
spinning motion in z-axis becomes 


^gg, X 


’ 3[-r sin ip Eo cos ( 0 + (o x t ) + (z - L c ) cos \j/ E J sin i cos w x t + r cos 6 

^gg, y 

= - n 2 

— 3[ — rsin ^ Eo cos(0 + cu z t) + (z - L c ) cos sin ip E „sin w x t + rsin0 

_ a gg.*_ 

spinning in z-axis 

- 3[ - r sin i > Eo cos (0 + io z t) + (z - L c ) cos 4/ E J cos \j/ Ea + (z - Q 


Substituting equations (4) with co x = o> z = 0 and 
equation (11) into equation (10), the noninertial frame 
expression of the gravity gradient acceleration with slew 
motion along the y-axis becomes 

3[ — r sin i// cos 0 + cos \ \>{z - L c )] sin^/r 4- rcosfl 
r sin 6 

- 3[ -r sin \p cos 6 4- cos \j/{z = L c )] cos \p + (z - L c ) 


fl gg ,y 

= —n 2 

_ a gg.zJ 

slew in y-axis 


( 12 ) 


(13) 


where ij/ E = ^ Eo 4- co y t. 

Substituting equations (4) with <o y = co z = 0 and 
equation (11) into equation (10), the noninertial frame 
expression of the gravity gradient acceleration with slew 
motion along the x-axis becomes 


a 

a 


gg* 

gg,y 


L gg- 2 . 


islew in x-axis 


3 A sin \j/ Ea 4- r cos 9 
3 A cos ip E sin co x t 4- r sin 0 
3 A cos \ j/ Eo cos w x t 4- (z — L c ) 


(14) 


where 


A = — r cos 0 sin ^ Eo 

4- cos \j/ E £r sin co x t 4- cos co x t(z — LJ] 

Thus, the gravity gradient acceleration at (r, 0, z) can be 
computed from that located at (x, y, z), shown in 
equations (12) to (14), from the following relation: 


fl gg,r 
^gg. o 



cos 0 
— sin 0 

sin 0 
cos 9 

o" 

0 

a gg,x 
fl gg. y 

J^gg.Z- 


0 

0 

1 

_ tt gg. 


3.5 Jitter accelerations 

In addition to the gravity gradient acceleration acting 
on the fluid element of on-orbit spacecraft fluid systems, 
there is another acceleration of gravity jitter that also 
exerts forces on fluid systems. Among all of the varieties 
of jitter accelerations listed, accelerations induced by slew 
motion of the spacecraft dominate over the forces 
activated on the spacecraft fluid systems. The spacecraft 
center of mass is located at (x c , y f , z c ) — (0, 0, L c ). 

A detailed expression of [R r , k d , & z ] slfW are shown in 
equation (9) of this paper. By putting [R r , R 0 , R z ] slew = 
[F r , F e , FJ slew , jitter acceleration can be expressed as 
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a th' 


fl gj.0 
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_°gj.z_ 
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00 

L__ 

olhers 




F r 
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[1 + 3 sin (2 7t/0] 
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F z 

slew 

1/.J 

others 



where / is the jitter frequency (Hz) imposed on the fluid 
systems of the spacecraft. 


4. Noninertial frame mathematical formulation of 
fundamental equations 

The dynamical behavior of fluids inside on-orbit 
spacecraft fluid systems are strongly affected by gravity 
gradient and gravity jitter accelerations. In order to 
accommodate the effect of gravity gradient acceleration 
on the on-orbit fluid motion, one has to consider a 
noninertial frame of the spacecraft-bound coordinate 
rather than an inertial frame coordinate usually adapted 
in ordinary fluid mechanics formulations. 

Consider a closed circular cylindrical dewar of radius, 
u, and height, L, which is partially filled with cryogenic 
liquid helium, and the rest of the ullage is filled with a 
helium vapor. The angular velocity of rotating cylinder 
is a>. Density and viscosity of liquid helium and helium 
vapor are p L , ^ L , and p v , respectively. 16,24,26 Let us 
use cylindrical coordinates (r, 0 , z), with corresponding 
velocity components (w, v, w), and corresponding residual 
gravity acceleration, such as gravity gradient compo- 
nents (flgg.r. flgg,* fl gg,z) and gravity jitter components 
(a gj a gj<0 , a gj z ). In the derivation of the governing 
equations, accelerations induced by the spinning motion 
of the spacecraft are included in the formulation. The 
governing equations for noninertial frame of spacecraft- 
bound coordinates can be shown as follows. 


4. 1 Continuity equation 

The continuity equation is 

id 1 dv dw 

(ru) + + — = 0 

r dr r dd dz 


( 17 ) 
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4.2 Momentum equations 
Radial direction : 

( du du v du v 

7 l + u 7 r + ;To~-r 

dp 


2 du 

+ W Jz 


= - ^ + p(a gj r + a gg r ) + praj 2 


f , u 2 
Circumferential direction : 


(18) 


with the characteristics of gravity gradient and gravity 
jitter accelerations as formulated in equations (1) to (16). 
Initial and boundary conditions shall be introduced to 
solve for the fluid motion in an on-orbit spacecraft fluid 
system in noninertial frame coordinates. 14,34-37 

Let the profile of the interface between gaseous and 
liquid fluids be given by 

niU r, M = 0 ( 23 ) 

The initial condition of the profile of the interface 
between gaseous and liquid fluids at t = t Q is assigned 
explicitly and is given by 


( dv dv v dv uv dv 

P\ (- U 1 1 h W — 

\dt dr r 36 r dz 


1 dp 
~ r 36 

+ H 

Axial direction: 


2 pco 2 u + p(a ihB + a gg # ) - prtb 1 
2 du > 




2 + r 2 se 


(19) 


r](t = t Q ,r, 6, z) = 0 (24) 

A set of boundary conditions has to be supplied for 
solving the equations. The initial interface profiles used 
in this study have been given explicitly through the 
steady-state computations made by Hung and Leslie 38 
and Hung and Shyu, 39,40 which were checked with the 
experiments carried out by Leslie. 15 These boundary 
conditions are as follows: 


( dw 3w v dw <3w 

V u h + w ----- 

dt dr r 36 dz 


( 20 ) 


where 


j 1 d ( d\ 1 d 2 d 2 
V = rJr\Fr) + ? 86 2 + 8? 


( 21 ) 


In these equations, 2 co z v and 2c o z u are the Coriolis 
acceleration, ra> 2 is the centrifugal acceleration, and 
rcb z is the angular acceleration induced by the spinning 
motion of the spacecraft. 

In the computation of fluid forces, moments, viscous 
stresses, and angular moments acting on the container 
wall of the spacecraft, one has to consider those forces 
and moments in the inertial frame rather than in the 
noninertial frame. To show an example, one has to 
transform those vectors from the noninertial frame to the 
inertia frame for the case of spinning motion in the z-axis. 
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cos c o z t 

— sin a ) z t 

o' 

F x 


F, 

= 

sin a) z t 

cos c o 2 t 

0 

Fy 

(22) 
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0 
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F t 



where primes denote vectors in the inertia frame 
while those parameters without the primes indicate 
vectors in the non-inertial frame. 


5. Initial and boundary conditions of spacecraft 
fluid system in the microgravity environment 

The governing equations of the fluid motion in on-orbit 
spacecraft fluid systems in a noninertial frame of 
spacecraft-bound coordinates have been illustrated in 
equations (17) to (22). These equations shall be combined 


(1) Along the container wall, the following three 
boundary conditions apply: 

(a) Interface between solid and liquid: No penetra- 
tion and no slip conditions assure that both the 
tangential and the normal components of the 
liquid velocity along the solid walls vanish. 

(b) Interface between solid and gaseous (vapor) 
fluid: Similar no penetration and no slip 
conditions as that for the interface between 
solid and liquid will apply. 

(c) At the location of solid-liquid-gaseous (vapor) 
three-phase interface: No penetration, but slip 
conditions apply. This will assure that normal 
components of liquid and vapor velocities 
along the solid wall vanish and allow a slipping 
flow of liquid and vapor fluids along the solid 
wall at the three-phase interface. The velocity 
of slipping flow at this location is governed by 
the adhesive forces between fluids (liquid and 
gas) and solid walls. Also, at the three-phase 
interface, a constant contact angle is present in 
which the behavior of wet or dry contacts are 
determined by molecular physics scale Cou- 
lomb interaction between the fluids (liquid and 
vapor) and the surface phenomena (material, 
polishing, and roughness) of solid walls. Mathe- 
matical computation of the interphase contact 
angle is very complicated, involves the molecular 
physics scope of microscale, and is not realistic 
to carry out such a computation in this sloshing 
dynamics modelling. In this study, experimental 
measurement of contact angle is employed in the 
mathematical modelling. 

(2) Along the interface between the liquid and gaseous 
fluids, the following two conditions apply: 

(a) Kinematic boundary condition : The liquid (or 
gaseous) surface moves with the liquid (or gas) 
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which implies 
Dr] 

— = 0 
Dt 

dr] dr] v dr] dt] 

dt dr r dO dz 

on rj(t = t a , r, 0, z) 


(25) 


(b) Interface stress conditions: At the liquid-vapor 
interface, the stress must be continuous. Based on 
Landu and Lifshitz, the stress across the 
liquid-vapor interface can be expressed as 

(Pg ~ Pi)* i ~ [(Tjj)o - 



where R j and R 2 are the radii of curvature. 

The expressions of radii of curvature R { and R 2 in 
cylindrical coordinates from differential geometry are 
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In equation (26) 
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is the viscous stress tensor, [i is the viscous coefficient 
of the first kind, f is the viscous coefficient of the second 
kind, P is the pressure, a is the surface tension of the 
liquid-vapor interface, is the unit vector normal to the 
interface, 6 {j is the Kronecker’s delta, and subscripts G 
and L denote conditions in the gaseous and liquids fluids, 
respectively, across the liquid-vapor interface. 

The fluid stresses across the liquid-vapor interface can 
be decomposed into a normal (n„ a unit vector) and 
tangential (f„ a unit vector) components to the interface. 
For the tangential component to the interface, one can 
take a dot product of a unit vector tangential to the 
interface, t h to equation (26), which yields 


C(Vi«j)] L = (28) 

since n^i = 0. 

For the normal component to the interface, one can 
also take a dot product of a unit vector normal to the 
interface, n h to equation (26), which leads to 



6. Methods of mathematical modelling simulation 

Detailed descriptions of the computational algorithm 
applicable to cryogenic fluid management under 
microgravity are also given in our earlier studies. 37,42 In 
this study, full-scale GP-B and AXAF spacecraft 
propellant dewar tanks with a radius of 68 cm and a 
height of 145 cm will be used in the numerical simulation 
(. Figure J). The propellant tank is 80% filled with 
cryogenic liquid helium and the rest of the ullage is filled 
with helium vapor. Temperature of the cryogenic helium 
is 1.3°K. In this study, the following data were used: 
liquid helium density = 0.145 g/cm 3 , helium vapor 
density = 0.00147 g/cm 3 , fluid pressure = 1.66 x 10 3 
dyne/cm 2 , surface tension at the interface between the 
liquid helium and the helium vapor = 0.346 dyne/cm, 
liquid helium viscosity coefficient = 1.12 x 10 _4 cm 2 /s 
(first viscosity coefficient only is considered because of 
the nature of incompressible flow which makes the effect 
of second viscosity vanishing in this case); and contact 
angle = 0°. The initial profiles of the liquid-vapor 
interface for the rotating dewar are determined from 
computations based on algorithms developed for the 
steady-state formulation of microgravity fluid manage- 
ment. 36,37 

A staggered grid for the velocity components is used 
in this computer program. The method was developed 
by Harlow and Welch 43 for their MAC (marker-and-cell) 
method for flows with a free surface. The finite difference 
method employed in this numerical study was the 
“Hybrid Scheme” developed by Spalding. 44 The 
formulation for this method is valid for any arbitrary 
interface location between the grid points and is not 
limited to middle point interfaces 45 A semi-implicit 
method 46 was used as the procedure for modelling the 
flow field. The time step is determined automatically 
based on the size of the grid space and the velocity of 



Figure 3. The GP-B module showing main elements of liquid 
helium dewar, probe, and baffle-boards. 
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Figure 4. Distribution of grid points in the (A) radial-axial plane 
and (B) radial-circumferential plane of the cylindrical coordinates 
for the GP-B dewar tank. 
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Figure 5. Distribution of grid points in the (A) radial-axial plane, 
and (B) radial-circumferential plane of the cylindrical coordinates 
for the AXAF-S dewar tank. 

the flow field. A detailed description of the computa- 
tional algorithm applicable to microgravity fluid 
management is illustrated in our earlier studies. 14 Figures 
4a and 4b show the distribution of grid points for the 
dewar tank with the probe for the GP-B dewar container 
in the radial-axial plane and radial-circumferential plane, 
respectively, in cylindrical coordinates. Figures 5a and 5b 
show the distribution of grid points for the dewar tank 
for the AXAF-S spacecraft in the radial-axial plane and 
radial-circumferential plane, respectively, in cylindrical 
coordinates. 

7. Spacecraft sloshing dynamics associated with 
spinning and/or slew motions 

By using the mathematical formulations illustrated in the 
previous sections, one can simulate numerically the 
spacecraft sloshing dynamics associated with spinning 
and/or slew motions depending upon the specific 
scientific missions assigned to the spacecraft. Examples 
are given to illustrate the sloshing dynamics associated 
with spinning motion for the GP-B spacecraft and the 


sloshing dynamics associated with slew motion for the 
AXAF-S spacecraft. 

7. 1 Sloshing dynamics associated with spinning motion 
for the GP-B spacecraft 

As described earlier, the dewar container of the GP-B 
is spinning with a rotation rate of 0.1 rpm during its 
normal operation. The noninertial frame equations, 
shown in equations (17)-(21), subject to initial and 
boundary conditions, shown in equations (23H29), shall 
be considered in our computation. As for gravity 
gradient acceleration, the noninertial frame expressions 
of equations (12) and (15) for spinning motion in the 
z-axis shall be adopted. For gravity jitter acceleration, 
all induced forces except forces due to slew motion shall 
be considered since the acceleration induced by the 
spinning motion has been included in the equations, 
shown in equations (18) and (19). Thus, gravity jitter 
acceleration, shown in equation (4), can be rewritten as 



sin \j/ Eo cos (0 + a) 2 t) 

— sin i/f En sin (0 + c o 2 t) {1 + \ (sin 27t/f))(30) 
cos 

Figure 6 shows the time variation of gravity gradient 
acceleration for a turnaround period of 1,200 s 
with container rotating speed (a> z ) of 0.1 rpm for the 
components along the (x, y, z) directions acted on the 
fluid mass located at (r, 0, z) = (40 cm, rc/4, 10 cm). Figure 
6 shows that the magnitude of gravity gradient 
acceleration is on the order of 10~ 7 g D ( g 0 = 981 cm/s 2 ). 
Figure 7 shows the time variation of the gravity jitter 
acceleration for a turn-around period of 1,200 s with a 
container rotating speed (cu z ) of 0.1 rpm and a jitter 
frequency (/) of 0.1 Hz for components along the (x, y, z) 
directions acted on the fluid mass everywhere in the 
container. Two ranges of background cavity, 10 6 and 
1O- 8 0 O f° r gravity jitter accelerations that correspond 
to accelerations higher and lower, respectively, than that 
of the gravity gradient acceleration are presented. 

The grid point distribution in the dewar tanks, shown 
in Figure 4 , is adopted for this example of the GP-B 
dewar. The present study is basically a simulation of an 
initial value problem. An improper assignment of initial 
profile of bubbles for the initial conditions will result in 
a wrong conclusion of simulation in this study. The 
initial conditions of the bubble shape shall be computed 
from the conditions of bubble configuration with 
minimum Gibbs free energy which is a function of surface 
tension, gravitational, and centrifugal accelerations 
together with the liquid-fill level and geometry of the 
container. The equilibrium shape of the liquid-vapor 
interface for a rotating dewar under a residual gravity 
environment below 1 O _ 6 0 O and a rotating speed of 
0.1 rpm is very much alike and is a doughnut 
configuration with a near circular kidney-shaped 
cross-section based on the computation of the numerical 
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Figure 6. Time variation of GP B spacecraft gravity gradient 
acceleration acting on fluid mass located at (r,6,z) = ( 40 cm, 
tt/ 4, 40 cm) for a turn-around period of 1 ,200 s with rotating speed 
of 0.1 rpm, along (A) the x-direction; (B) the y-direction; (C) the 
z-direction. w = 0.1 rpm; r = 1 ,200 s; g 0 = 9.81 m/s 2 . 


algorithm developed in our earlier studies, 34,35 Figure 8a 
shows the initial profile of the liquid-vapor interface in 
the r-z plane at 9 = 0° and 180°. Figure 8b shows the 
initial profile of the liquid-vapor interface in the r-z 
plane at 6 = 90° and 270°. Figure 8c shows the initial 
profile of the liquid-vapor interface in the r-6 plane at 
a height z = 108 cm. Figure 8d shows the initial profile 
of the three-dimensional liquid-vapor interface. By using 
initial conditions of the bubble profiles described in 
Figure 8 , one can proceed with the simulation of bubble 
sloshing dynamics activated by orbital accelerations. 

The sloshing dynamics of the GP-B spacecraft driven 
by the combined effects of gravity gradient and a 10“ 8 g Q 
background jitter acceleration associated with spinning 
motion in the z-axis has been studied. In this case, the 
combined forces are dominated by the gravity gradient 
acceleration. Figure 9 shows the time sequence evolution 
of the three-dimensional dynamical behavior of the 
interface oscillations driven by these combined accelera- 
tions. The figure shows the time sequence of the bubble 
oscillations due to the GP-B sloshing dynamics at times 


Hung and H. L Pan 
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Figure 7. Time variation of GP-B spacecraft gravity gradient 
acceleration acting on the fluid system under background gravity 
of 10“ 8 and 10' 6 g ot rotating speed of 0.1 rpm, turn- around period 
of 1,200 s and gravity jitter frequency of 0.1 Hz, along (A) the 
x-direction; (B) the y-direction; and (C) the z-direction. 


t = 191, 354, 380, 431, 503, 603, 825, 980, 995, 1,050, 
1,080, and 1,200 s. It clearly shows that there are a series 
of asymmetric oscillations excited along the liquid-vapor 
interface driven by the asymmetric gravity gradient 
dominated acceleration associated with spinning motion 
along the z-axis. 

With reference to equation (13), the gravity gradient 
acceleration exerted on the spacecraft is equivalent to the 
combination of time-dependent force with tidal motion 
and torsional moment acting on the spacecraft fluid mass 
when it is orbiting around the Earth. 6,7 Figure 9 shows 
the simulation of time animation of how the bubble 
deforms in response to the time-dependent force with 
tidal motion and torsional moment characterized by the 
gravity gradient dominated accelerations activated on 
the rotating dewar. 

Time fluctuations of the locations of bubble mass 
centers of the rotating container due to slosh waves 
excited by gravity gradient dominated acceleration based 
on Figure 9 have been computed. The values of bubble 
mass center fluctuations are (Ax c , Ay c , Az c ) — (4.54, 
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Figure 8. Initial profiles of liquid-vapor interface for the Gravity 
Probe- B module of a rotating dewar under background gravity of 
10“ 7 g Q , rotating speed of 0.1 rpm, direction of background gravity 
at \j/£ = 0°, and gravity jitter frequency of 0 1 Hz. In (A) the r-z plane 
at 8 = 0‘ and 180 , (B) the r-0 plane at 0 = 90 r and 270 , (C) the 
r-0 plane at z = 108 cm, and (D) the three-dimensional liquid-va- 
por interface profile. 


2.83, 2.32) cm. It shows that Ax c > Ay c > A z c for 
bubble mass center fluctuations driven by gravity 
gradient dominated acceleration. 

The combined effects of gravity gradient and a 1 0 6 g Q 
background jitter acceleration associated with spinning 
motion have also been studied. In this case, the combined 
forces are dominated by the jitter acceleration. 
Figure 10 shows the evolution of the three-dimensional 
dynamical behavior of the sloshing dynamics governed 
interface oscillations driven by the gravity jitter 
dominated acceleration of 0.1 Hz, i.e., a low jitter 
frequency. It clearly shows that there are a series of 
asymmetric oscillations excited along the surface of the 
liquid-vapor interface driven by asymmetric gravity jitter 
dominated acceleration. 

With reference to equation (17), gravity jitter 
acceleration exerted on the spacecraft is equivalent to 
time-dependent oscillatory forces that push the bubble 
in the combined directions of down-and-up and 
leftward-and-rightward as the bubble is rotating with 
respect to the spacecraft rotating axis. Figure 10 shows 
the simulation of how the bubble deforms in response to 
the time-dependent up-and-down and leftward-and- 
rightward forces characterized by the jitter dominated 
accelerations activated on the rotating dewar. 

Time fluctuations of the locations of bubble mass 
centers of the rotating container due to slosh waves 



Figure 9. Time sequence evolution of the GP-B spacecraft 
three-dimensional liquid-vapor interface oscillations for a rotating 
dewar driven by combined gravity gradient and jitter accelerations 
with background gravity of 10 -a g o . The rotating speed of the 
dewar is 0.1 rpm, and the period of gravity direction turn-around 
time is 1 ,200 s. 


excited by asymmetric jitter-dominated acceleration 
based on Figure 10 have been computed. The values of 
bubble mass center fluctuations are (Ax c , Ay c , Az c ) = 
(17.28, 12.2, 33.83) cm. A comparison of Figures 9 and 10 
for bubble mass center fluctuations between that driven 
by gravity gradient dominated and that given by jitter 
dominated accelerations provides the following conclu- 
sions. (a) Dynamics of bubble (liquid-vapor interface) 
driven by gravity gradient dominated acceleration 
produces bubble mass center fluctuations of Ax c > 
Ay c > Az c while up and down oscillations of bubble 
driven by jitter dominated acceleration induces bubble 
mass center fluctuations of Az c > Ax c > Ay c . (b) A 
comparison of the values of bubble mass center 
fluctuations show that fluctuations of bubble mass center 
induced by jitter dominated acceleration is greater than 
that induced by gravity gradient dominated acceleration, 
(c) Fluctuations of both x c and y c start from zero while 
that of the z c starts from a non-zero value at the middle 
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Figure 10. Time sequence evolution of the GP-B spacecraft 
three-dimensional liquid-vapor interface oscillations driven by 
combined gravity gradient and jitter accelerations of 0.1 Hz 
frequency, background gravity of 10~ 6 g 0 with the period of gravity 
direction turn-around time = 1,200 s, and rotating speed of 
0.1 rpm. 


point of the height of partially filled liquid container, (d). 
Fluctuation magnitudes of x c , y c , and z c associated with 
higher jitter frequency are lower than those of the 
frequencies associated with lower jitter frequency for 
bubble mass center fluctuations driven by jitter 
dominated acceleration. 

7.2 Sloshing dynamics associated with slew motion for 
the AXAF-S spacecraft 

Assuming that the slew motion is along the /'-axis 
(Figure 2), the gravity gradient acceleration associated 
with slew motion can be computed from the noninertial 
frame expressions of equations (8) and (10). It is assumed 
that the slew motion operates at 90° in 600 s. 

In this example, the spacecraft sloshing dynamics 
driven by gravity gradient acceleration associated with 
slew motion along the /'-axis, shown in Figure 2, have 
been investigated. As the orbital period of the AXAF-S 
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Figure 11 . Time variation of AXAF-S spacecraft gravity gradient 
acceleration acting on fluid mass located at {r, 6, z) — (12 cm, 
n/2, 3 cm) for 90° slew motion in 1 0 minutes along the /'-axis and 
orbital period of 97.6 minutes along (A) the x-direction; (B) the 
y-direction; (C) the z-direction. g Q = 9.81 m/s 2 . 


spacecraft is 97.6 minutes, and the period of slew motion 
along the /'-axis is 600 s. The component of gravity 
gradient acceleration along the (x, y , z) directions on a 
fluid mass located at (r, 0, z) = (12 cm, rc/2, 3 cm) is 
shown in Figure 11. This figure shows that the magnitude 
of the gravity gradient acceleration is on the order of 
10“ 7 0 O for the AXAF-S dewar on its operation orbit. 
The distance from the spacecraft mass center to the 
bottom of the dewar, L c , shown in Figure 2, is 257.8 cm. 
The grid point distribution of the dewar tanks, shown in 
Figure 5, is adopted for the AXAF-S dewar. 

As indicated earlier, the initial shape of the bubble 
shall be given with the consideration of surface tension 
and gravitational acceleration, together with the liquid 
fuel level and the geometry of the container based on the 
computation of the minimum Gibbs free energy of the 
bubble. The equilibrium shape of the liquid-vapor 
interface for a dewar with 70% liquid-filled level under 
a residual gravity environment below 10“ 7 g a is a sphere. 
Figure 12a shows the initial shape of the interface in the 
r-z plane at 0 — 0° and 180°. Figure 12b shows the initial 
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(A) r-z Diane at Q * 0* and 180* (B) r-z Diane at 0* 90* and 180* 
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Figure 12. Initial profiles of liquid-vapor interface for AXAF-S 
spacecraft of dewar tank under background gravity of 1 0 7 g QI and 
direction of background gravity at i \j/ t = 0 In (A) the r-z plane at 
9 = 0° and 180°, (B) the r-6 plane at 0= 90 and 270°, (C) the r-9 
plane at z = 95.9 cm, and (D) the three-dimensional liquid-vapor 
interface profile. 


profile of the liquid-vapor interface in the r-z plane at 
0 — 90° and 270°. Figure 12c shows the initial profile of 
the liquid-vapor interface in the r-0 plane at height 
z = 95.9 cm. Figure I2d shows the initial profile of 
three-dimensional liquid-vapor interface. Time-depen- 
dent simulation can be performed based on the initial 
conditions with the bubble profiles specified in Figure 12. 

The sloshing dynamics of the AXAF-S spacecraft 
driven by the gravity gradient acceleration associated 
with slew motion along the /'-axis, shown in Figure 2, 
have been investigated. Figure 13 shows the time 
evolution of the three-dimensional dynamical behavior 
of the interface oscillations driven by gravity gradient 
acceleration associated with slew motion. 

With reference to the characteristics of gravity 
gradient acceleration associated with slew motion, shown 
in Figure 77, liquid is pushed toward the negative 
x-direction and the negative z-direction with gradually 
decreasing values initially and increasing values later but 
still in the negative direction. This makes the bubble to 
be pushed toward the positive x-direction and positive 
z-direction with decreasing and increasing accelerations 
correctively. The uneven and imbalance flow velocities 
induced by gravity gradient acceleration associated with 
slew motion toward the positive x- and positive 
z-directions create a similar uneven and imbalanced 
pressure distribution on the bubble which have resulted 



Figure 13. Time sequence evolution of the AXAF-S spacecraft 
three-dimensional liquid-vapor interface oscillations for a dewar 
driven by gravity gradient acceleration associated with slew motion 
in the /’-axis. An orbital period of 97.6 minutes and 90 slew motion 
in 10 minutes are applied to the spacecraft operation. 


in a deformed irregular concave- and convex-shaped 
oscillating bubble shown in Figure 13. 

Time fluctuations of the locations of bubble mass 
centers of the fluids inside the dewar container due to 
sloshing dynamics driven by gravity gradient accelera- 
tion associated with slew motion based on Figure 13 have 
been computed. The values of bubble mass center 
fluctuations are (Ax c , Ay c , Az c ) = (1 1.8, 1.25, 18.4) cm. It 
shows Az c > Ax c > Ay c for bubble mass center fluctua- 
tions driven by gravity acceleration associated with slew 
motion. 

The trend of the bubble mass center fluctuations, 
Az c > Ax c > A y c , exactly reflects the values of the major 
driving forces of gravity gradient acceleration associated 
with slew motion in which it shows g z > g x > g y as that 
shown in Figure 77. 

For a slew motion that operates with a range of 9CF 
in 600 s, the component of jitter acceleration, based on 
equations (6), (8), (9), and ( 1 7), along the (x, >\ z) directions 
acted on the fluid mass located at (r, 0, z) = (12 cm. 
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7t/2, 3 cm) is shown in Figure 14. This figure shows that 
the magnitude of gravity gradient acceleration is on the 
order of 10~ 5 g Q . 

Figure 15 shows the time evolution of the three- 
dimensional behavior of the interface oscillations driven 
by jitter acceleration associated with slew motion. In this 
figure, a time evolution of liquid-vapor interface profiles 
is shown at times t = 60.4, 258, 323, 354, 379, 403, 430, 
528, 572, 628, 663, and 800 s. With reference to the 
characteristics of jitter acceleration associated with slew 
motion, shown in Figure 14 , liquid is pushed toward the 
positive x- first and then to the negative x-directions, in 
combination with the negative z-direction. This makes 
the bubble to be pushed toward the negative x- first and 
then to the positive x-direction, together with the positive 
z-direction. These fairly complicated time-dependent 
variations of the dynamical forces exerted on the bubble 
create a series of bubble deformations as shown in 
Figure 15. 

Time evolution of the locations of bubble mass centers 
of fluids inside the dewar container due to sloshing 
dynamics driven by jitter acceleration associated with 
slew motion based on Figure 15 have been computed. 




Figure 14. Time variation of AXAF-S spacecraft jitter acceleration 
associated with slew motion acting on fluid mass located at 
( r ; 6 , z) = (1 2 cm, n/2, 3 cm) for 90° slew motion in 600 s along the 
/'-axis and orbital period of 97.6 minutes along (A) the x-direction; 
(B) the /- direction; and (C) the z-direction. 



Figure 15. Time sequence evolution of AXAF-S spacecraft 
three-dimensional liquid-vapor interface oscillations for a dewar 
driven by jitter acceleration associated with slew motion in the 
/'-axis. 90° slew motion in 600 s is applied to the spacecraft 
operation. 


The values of bubble mass center fluctuations are 
(Ax c , Ay c , Az c ) = (28.9, 0.44, 30.2) cm. It shows Az c > 
Ax c > A y c for bubble mass center fluctuations driven by 
jitter acceleration associated with slew motion. Behavior 
of bubble mass center fluctuations, shown in Figure 75, 
are the exact reflection of the behavior of jitter 
accelerations, shown in Figure 14. 

Comparison between Figures 14 and 15 illustrate 
some peculiar behavior of helium fluids with temperature 
below A-point (2.17°K) in which helium demonstrates a 
number of remarkable properties of superfluidity, such 
as extremely low viscous and surface tension coefficients 
reacted to the disturbances driven by jitter acceleration 
associated with slew motion. It can be concluded as 
follows: (a) jitter acceleration associated with slew 
motion started at time t = 0. However bubble mass 
center fluctuations did not start to react to the driven 
force until t = 80 s for x c and t — 160 s for z c . In other 
words, the bubble configuration was perfectly spherical- 
shaped at 80 s after jitter acceleration associated with 
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slew motion was applied, (b) The x-component of jitter 
acceleration associated with slew motion was applied to 
the fluid element with a positive value from t — 0 to 
300 s; a negative value from t = 300-600 s; and zero 
value after t = 600 s. However, westward movement of 
the bubble continued from t = 80-450 s and then 
switched to eastward movement to t = 800 s, which was 
200 s after the acceleration vanished. Obviously there is 
a phase shift between the action of force and the reaction 
of motion, (c) The z-component jitter acceleration 
associated with slew motion was applied to the fluid 
element with a negative value from r = 80-500 s and the 
zero value after t = 500 s. However, the northward 
movement of the bubble started at t = 160 s, bounced 
back from the wall at t = 530 s, and then southward 
movement continued to t = 800 s, which w r as 300 s after 
the acceleration vanished. Obviously the motion 
continued for a long period of time even after the applied 
force vanished due to extremely low viscosity of helium 
II fluids, (d) An intensive oscillation of bubble with a 
deformation of irregular concave and convex shaped 
continued and sustained for several hundred seconds 
after the applied force vanished due to an extremely low 
surface tension coefficient for helium II fluids. 


8. Discussion and conclusion 

For different scientific missions, some experimental 
spacecrafts have to operate with various kinds of slew 
motion to perform their scientific experiments. In this 
study, both the GP-B and AXAF-S spacecrafts have been 
used as examples to show the mathematical modelling 
of bubble sloshing dynamics for cryogenic liquid helium 
in orbital spacecraft dewar containers. The instability of 
the liquid-vapor interface of the fluid systems can be 
induced by the presence of gravity gradient and/or jitter 
accelerations associated with spinning and/or slew 
motions of the spacecraft. These instabilities produced 
by sloshing dynamics can cause various problems in 
spacecraft control systems. Sometimes these sloshing 
dynamics can even deteriorate the quality of the normal 
operation of the spacecraft. 

In this paper, a general mathematical model of 
sloshing dynamics of the dewar container partially filled 
with liquid driven by both a gravity gradient and jitter 
accelerations associated with spinning and/or slew 

motions were investigated. Two examples of the 
numerical simulation of sloshing dynamics driven by 
orbital accelerations, one with the GP-B for spinning 
motion and the other with the AXAF-S for slew motion, 
have been carried out. These examples demonstrate the 
variety of sloshing dynamics caused by the orbital 
accelerations associated with spinning and/or slew 

motions. 
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